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FOREWORD 

The research results reported here were prepared for the National 
Aeronautics and Space Administration, Johnson Space Center (NASA, 
JSC) under a grant to the Environmental Research Institute of Michi- 
gan (ERIM), Radar and Optics Division. The technical monitor for 
NASA, JSC was Dr. M. B. Duke. Dr. P.L. Jackson* was Principal In- 
vestigator during the first 12 months of the program; Mr. R. F. Raw- 
son was designated Principal Investigator for the completion of the 
program. R.E. Hamilton and C.L. Liskow provided analysis support 
and A.R. Dias provided support in using the ERIM digital processing 
facilities. 


*Dr. Jackson is presently with the University of Michigan, De- 
partment of Geology, Ann Arbor, Michigan. 
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INTRODUCTION AND OBJECTIVES 


1.1 INTRODUCTION 

Data obtained with synthetic aperture radar (SAR) is processed 
to obtain imagery with intensity proportional to the reflectivity 
distribution of the scene imaged. During the past decade, and more, 
a large number of investigators have studied ways to model the ob- 
served reflectivity distribution in terms of physical surface pro- 
perties for a variety of scattering surfaces. In order to validate 
models, measurements of the reflectivity of specific surfaces must 
be compared to the theoretical values of reflectivity as predicted 
by the models using ground measurements. In general practice, how- 
ever, the extensive ground measurements of specific imaged sites 
needed to make these validations have been both difficult and expen- 
sive to obtain. Work on correlating the reflectivity distributions 
with physical surface parameters through scattering models has there- 
fore been extremely limited. 

The work described in this report has made use of available SAR 
data obtained at the SP Mountain lava flow. This particular site 
has been the subject of numerous investigations which have provided 
a wealth of scattering model information based on extensive ground 
measurements. The research reported here has taken advantage of both 
the SAR data collected by ERIM with the ERIM X-L SAR system, and the 
available model information, to correlate reflectivity distributions 
obtained from SAR data with the ground measurements and models de- 
scribing the imaged area. The correlations were accomplished for 3 
cm wavelength data (X-band) and 23 cm wavelength data (L-band). Both 
like (horizontal) and cross polarization data sets were utilized. 
Various image manipulations are required in order to fully utilize 
the SAR data for scattering studies. The availability of the X-L 
SAR data in particular has allowed the opportunity to demonstrate 
several sophisticated digital processing techniques. 
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1.2 OBJECTIVES 

The research effort discussed in this report was carried out for 
the general purpose of analyzing existing X-L radar data of SP Moun- 
tain and its associated lava flow in such a way as to correlate the 
radar return with terrain surface roughness. In particular, the ob- 
jectives of this analysis were to (1) analyze all four channels in 
terms of geological applications, (2) determine if the vegetation 
masks the geology, (3) apply Multispectral Scanner (MSS) enhancement 
techniques to the radar imagery, and (4) combine radar and MSS chan- 
nels together. 

The following chapters of this report discuss the steps taken 
toward meeting these objectives- Chapter 2 describes the data 
processing carried out to prepare the radar imagery for analysis and 
interpretation. Chapter 3 presents the results of manual interpre- 
tation of the imagery and discusses the digital image processing 
techniques applied to aid geologic interpretation. Chapter 4, the 
final chapter, presents the conclusions of the analysis and recommen- 
dations for further research in this area. 
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PREPARATION OF RADAR IMAGERY 

This analysis made use of radar image data collected on 20 Octo- 
ber 1979. During Pass 3 on that date, the SP Mountain area in Ari- 
zona was imaged by means of ERIM's X-L four-channel imaging radar 
mounted in a Canadian CV-580 aircraft. The flight altitude (i.e., 
height above ground) was approximately 15,000 ft and the antenna de- 
pression angle was 36°. The system used two transmitter channels, 
X-band and L-band, both at horizontal polarization. There were four 
receiver channels, two X-band and two L-band; for each wavelength, 
one channel received parallel-polarized signals and the other re- 
ceived cross-polarized signals*. 

Signal data from the X-L SAR is recorded on film on-board the 
aircraft. Imagery is obtained from the signal data using a matched 
filter type processor in the ERIM optical processing facility. Two 
image formats can be obtained for analysis: (1) a film transparency 
and (2) digital data on a CCT. In this investigation, film trans- 
parencies of the imagery were utilized only to identify the test area 
to be used for further analysis. All subsequent data reduction was 
accomplished using the images in a CCT format. 

This section describes, briefly, the steps in preparing the radar 
imagery for analysis. The preparations include the following steps: 
(1) digitization of the optically-processed imagery, (2) image gen- 
eration and annotation, (3) image registration, (4) slant range-to- 
ground range coordinate transformation, (5) scale corrections, and 
(6) image intensity corrections. These various steps are described 
in this section. 



*R.F. Rawson, F. Smith, and R. Larson, "The ERIM Simultaneous 
X- and L-Band Dual Polarization Radar", The Record of the IEEE 1975 
International Radar Conference, 21-23 April 1975. 
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2.1 DIGITIZATION 

The signal film which resulted from the data-col lection mission 
was used as the input to ERIM's Hybrid Image Processing Facility 
(HIPF). This facility converted the analog image data to digital 
form and recorded the digital data on computer-compatible tape (CCT). 
Parameters of the digitization process are discussed in Appendix A. 
The resulting tapes, designated IPL 584, 585, 586, and 587, served 
as inputs to the data-analysis task discussed in this report. 

2.2 IMAGERY GENERATION 

The four (raw data) CCTs produced by HIPF were used by ERIM's 
Earth Resources Data Center (ERDC) to generate the uncorrected imag- 
ery shown in Figures la, 2a, 3a, and 4a. These radar images were 
all printed with north down in order to have the computer-printed 
identification oriented for easy reading. 

Early in this program, it was decided that all imagery should be 
carefully labelled in order to avoid confusion. Therefore, a pro- 
cedure was developed for placing pertinent information on each trans- 
parency along with the image. Using the CCT as its input, ERDC 
generated an image transparency upon which appears the header infor- 
mation plus certain letter codes which indicate the correction 
process applied to the data. All imagery presented in this report 
is so labelled. 

2.3 IMAGE REGISTRATION 

In order to compare the various radar-channel images with one 
another using a digital computer, it is necessary to register all 
four images. This registration requires a geometric reference, such 
as a topographic map, a satellite photo, or an appropriate aerial 
photo. 


4 
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(a) Uncorrected 


(b) SR-to-GR Corrected 


FIGURE 1. ERDC IMAGERY, X-BAND, HH POLARIZATION. 





(a) Uncorrected (b) SR-to-GR Corrected 


FIGURE 2. ERDC IMAGERY, X-BAND, HV POLARIZATION. 
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(a) Uncorrected (b) SR-to-GR Corrected 


FIGURE 3. ERDC IMAGERY, L-BAND, HH POLARIZATION. 
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(a) Uncorrected 


(b) SR-to-GR Corrected 


FIGURE 4. ERDC IMAGERY, L-BAND, HV POLARIZATION. 
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For this program, registration was accomplished in the following 
manner. First, since the original digital radar image was linear in 
slant range, a conversion to ground-range coordinates was accom- 
plished. Second, the range scale was adjusted with respect to the 
azimuth scale to obtain unity aspect ratio. These two steps were 
carried out in one operation using ERIM's Advanced Radar Imaging and 
Exploitation System (ARIES) digital computer facility and are de- 
scribed in the next two sections. Details of the digital computer 
work are presented in Appendix B. 

2.3.1 SLANT-RANGE-TO-GROUND-RANGE CORRECTION 

To accomplish the slant-range-to-ground-range (SR to GR) cor- 
rection of the digital image data, an existing ERIM computer program 
was employed. This program employs a digital reconstruction filter 
to generate a new unequally-spaced set of data samples from the ori- 
ginal samples which are equally spaced in slant range; Figure 5 
illustrates the use of this filter. In order to effect the desired 
correction, the new samples must be unequally spaced in a very par- 
ticular way, namely, such that they correspond to equal spacings in 
ground range. This is accomplished by considering the geometrical 
relationship between the slant range and ground range coordinate 
systems, as shown in Figure 6. From this figure, the following equa- 
tion can be derived; 


- h . (1) 

where d^^ is a slant-range distance, d^^ is the corresponding 
ground-range distance, and h is the aircraft's altitude above the 
terrain; these terms are clarified in Figure 6. Using this equation, 
equal ground range spacings ad^^ will result in the proper unequal 
slant range spacings Ad^^^ 


INPUT SAMPLES, s(t 



RECONSTRUCTION OF SLANT RANGE DATA. 







FIGURE 6. RELATIONSHIP BETWEEN SLANT RANGE AND GROUND RANGE COORDINATE SYSTEMS. 
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The imagery resulting from this SR-to-GR correction process is 
shown in Figures lb, 2b, 3b, and 4b. Note that the code letter "G" 
(for "ground") has been appended to the header information to signify 
that this correction has been made. 

2.3.2 RANGE-AZIMUTH SCALE CORRECTION 

The original intent of this correction was that the range and 
azimuth scales on all four images be made equal. Unfortunately, this 
approach turned out to be very time-consuming and costly; it requires 
about eighteen hours of computer time per tape. An alternate 
approach was therefore selected which requires about one hour of com- 
puter time per tape; this approach corrects the ratio of the range 
and azimuth scales to unity, but does not make all four images exact- 
ly the same size. Further information about the two alternatives is 
given in Appendix B. 

In order to determine the required correction factors, measure- 
ments were made on the SR-to-GR-corrected images in Figures lb, 2b, 
3b, and 4b. The procedure for determining the range-to-azimuth scale 
ratio was to first measure the range distance between two identifi- 
able points spaced far apart in the radar image and to measure the 
corresponding distance on a satellite color photograph.* By di- 
viding the former distance by the latter, a range quotient was cal- 
culated. Then, the same types of measurements were made in the azi- 
muth direction and the corresponding azimuth quotient was formed. 
Finally, the range quotient was divided by the azimuth quotient to 
give the desired range-to-azimuth scale ratio. By simply inverting 
this ratio, the correction factor is obtained. Table 1 lists these 
factors for all four channels. It can be shown that these correction 


*This reference was chosen because it was easier to correlate 
the boundaries of the cone and lava flow between it and the radar 
image (than it was using the topographic map). 
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TABLE 1 

SP MOUNTAIN RANGE-AZIMUTH SCALE CORRECTION FACTORS 


Channel 

IPL 

File 

No. 

Pixel Spacing (m) 
Azimuth Range 

Scale 

Mul tipi ication 
Factor 

Azimuth Range 

Xhh 

584-1 

1.152 

1 

1.152 

>‘hv 

585-1 

1 R 1-5 

TTTT? 

1 

1.112 

'"HH 

587-1 

15 ^ 

1.308 

1 

1.308 

'-HV 

586-1 

1 ^'5 

T3U6 

1 

1.306 
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factors are within 7 percent of those indicated for the original 
ideal approach. The decision was made that it would not be cost- 
effective to use eighteen times as much computer time to eliminate 
an error of this magnitude. 

Images which have been both SR-to-GR-corrected and range-azimuth- 
scale-corrected were produced. Measurements were made on these 
images in order to determine how closely the desired unity range- 
azimuth scale ratio was achieved. The results are shown in Table 2. 
It can be seen that the errors in scale ratio range from 2.1 to 5.9 
percent. 

2.4 INTENSITY CORRECTION 

The intensity of the radar imagery has variations within the 
image due to (1) radar antenna gain variation, (2) radar range varia- 
tion, and (3) radar display irregularities. The radar image inten- 
sities can be corrected for antenna gain variation by increasing the 
intensity of image components at each depression angle to the value 
they would have with peak antenna gain. This is done by multiplying 
each value by the square of the ratio of peak gain to actual gain. 
All image components away from the peak gain angle will be increased 
in intensity. 

The radar image intensity variation due to radar range variation 
can be corrected by a factor 

-3 

PR 1 



where R^ is the range to the center of the scene. The received 
echo from identical targets at different ranges varies as R”^, but 
the compression factor attained by the synthetic aperture radar 
processor increases in proportion to R. These two variations with 
range result in the R relationship. Since the center range of 
the display was chosen as the reference range, the near-range 


14 


RADAR AND OPTICS DIVISION 



TABLE 2 

SCALE RATIOS ACHIEVED BY CORRECTION PROCESS 


IPL 

Number 

Radar 

Channel 

Scale 

Ratio 

Achieved 

584 

^HH 

1.027 

585 

^HV 

1.021 

587 


1.043 

586 

•-HV 

1.059 
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intensities will be decreased and the far-range intensities will be 
increased. 


Correction for radar signal recorder irregularities is accom- 
plished by use of a calibration curve generated from recordings of 
standard signals and measurements of the image intensities of these 
standard signals. All image component intensities are multiplied by 
the ratio of peak response to the response at the range of interest. 
This brings all image intensities up to the same efficiency across 
the radar display. 

Correction factors for these three intensity variations were 
calculated at 1000 ft intervals in slant range across the swath. 
These are summarized in Tables 3, 4, 5, and 6; the pertinent geometry 
is illustrated in Figure 7. Accurate interpolation at the far edge 
of the swath required that correction factors be provided out to 
40,000 ft. 

With these three intensity corrections, each intensity value can 
be converted to equivalent radar cross section by use of a cali- 
bration chart. This equivalent radar cross section could be useful 
for point reflectors, such as spheres or corner reflectors. However, 
extended horizontal surfaces will appear weaker at lower depression 
angles because the reflected power is inversely proportional to the 
cosine of the depression angle. This is due to the fact that equal 
sections of horizontal surfaces intercept less of the illuminating 
radar signal at smaller depression angles. Thus, the reflected sig- 
nal from extended targets is less at lower depression angles even 
though the surface reflectivity is constant. 

If different channels of radar imagery are to be compared by 
ratioing or subtraction, data corrected as described above would 
provide good results only for point targets. If, however, it is 
desired to compare imagery of horizontal surfaces, the intensities 
should be corrected for illumination differences. That is, they 
should be multiplied by cos e. 


• V: 
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TABLE 3 

INTENSITY CORRECTION INFORMATION 


Ground 

Range 

(ft) 

Percent 

of 

Display 

Width 

SI ant 
Range 
(ft) 

Depression 

Angle 

(Degrees) 

.,2 '“='0 [ S ;] 

(dB) (dB) 

! 

^eff 

(dB) 

Total 

Error 

(dB) 

6,213 

0 

16,236* 

67.5 


-5.75 



8,002 

4 

17,000 

61.92 


+5.65 



9,951 

10 

18,000 

56.44 


+4.40 

-0.2 


11,661 

15 

19,000 

52.14 


+3.70 

-0.2 


13,229 

21 

20,000 

48.59 

-14.5 

+3.03 

-0.2 

-11.67 

14,698 

26 

21,000 

45.58 

-11.0 

+2.39 

-0.2 

-8.81 

16,092 

32 

22,000 

42.99 

-9.6 

+1.79 

-0.1 

-6.91 

17,434 

38 

23,000 

40.71 

-6.7 

+1.21 

-0.0 

-5.49 

18,736 

43 

24,000 

38.68 

-4.7 

+0.65 

-0.1 

-4.15 

20,000 

49 

25,000 

36.87 

-2.8 

+0.12 

-0.6 

-3.28 

20,294 

50 

25,236 

36.47 

-2.5 

0.00 

-0.8 

-3.30 

21,232 

54 

26,000 

35.23 

-1.7 

-0.39 

-1.1 

-3.19 

22,450 

60 

27,000 

33.75 

-0.8 

-0.88 

-1.7 

-3.38 

23,644 

65 

28,000 

32.39 

-0.3 

-1.35 

-2.4 

-4.05 

24,821 

71 

29,000 

31.14 

-0.08 

-1.81 

-2.2 

-4.09 

25,981 

76 

30,000 

30.00 

0.0 

-2.25 

-2.1 

-4.35 

27,129 

82 

31,000 

28.94 

-0.05 

-2.68 

-2.0 

-4.73 

28,261 

86 

32,000 

27.95 

-0.23 

-3.09 



29,393 

93 

33,000 

27.04 

-0.53 

-3.49 



30,522 

99 

34,000 

26.12 

-0.95 

-3.88 



30,716 

100 

34,236 

25.98 

-1.0 

-3.97 



31,623 


35,000 

25.38 

-1.1 

-4.25 


-6.2 

32,726 


36,000 

24.62 

-1.8 

-4.63 


-7.2 

33,823 


37,000 

23.92 

-2.0 

-4.99 


-7.8 

34,914 


38,000 

23.25 

-2.5 

-5.33 


-8.6 

36,000 


39,000 

22.62 

-2.8 

-5.67 


-9.3 

37,081 


40,000 

22.02 

-3.2 

-6.00 


-10.0 


*Near edge. 
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TABLE 4 

INTENSITY CORRECTION INFORMATION 


Percent 

of 

Display 

Width 

SI ant 
Range 
(ft) 

Depression 

Angle 

(Degrees) 

(d8) 

^°9l0 [y 

(dB) 

"eff 

(dB) 

Total 

Error 

(dB) 

0 

16,236 

67.50 

-34.0 

5.75 

-5.0 

-33.3 

4 

17,000 

61.92 

-34.0 

5.65 

-5.0 

-33.4 

10 

18,000 

56.44 

-23.8 

4.40 

-5.5 

-24.9 

15 

19,000 

52.14 

-19.6 

3.70 

-5.6 

-24.5 

21 

20,000 

48.59 

-16.8 

3.03 

-5.9 

-19.7 

26 

21,000 

45.58 

-14.4 

2.39 

-6.0 

-18.0 

32 

22,000 

42.99 

-11.9 

1.79 

-5.8 

-15.9 

38 

23,000 

40.71 

-8.9 

1.21 

-5.4 

-13.1 

43 

24,000 

38.68 

-5.9 

0.65 

-5.1 

-10.4 

49 

25,000 

36.87 

-3.4 

0.12 

-5.0 

-8.3 

50 

25,236 

36.47 

-2.9 

0.00 

-5.0 

-7.9 

54 

26,000 

35.23 

-2.1 

-0.39 

-5.0 

-7.5 

60 

27,000 

37.75 

-0.9 

-0.88 

-4.8 

-6.6 

65 

28,000 

32.39 

-0.3 

-1.35 

-3.9 

-2.9 

71 

29,000 

31.14 

0.0 

-1.81 

-3.0 

-4.8 

76 

30,000 

30.00 

0.0 

-2.25 

-2.0 

-4.3 

82 

31,000 

28.94 

-0.0 

-2.68 

0.0 

-2.7 

86 

32,000 

27.95 

-0.2 

-3.09 

0.0 

-3.3 

93 

33,000 

27.04 

-0.7 

-3.49 

0.0 

-4.2 

99 

34,000 

26.12 

-0.9 

-3.88 

0.0 

-4.8 

100 

34,236 

25.98 

-1.1 

-3.97 

0.0 

-5.1 


35,000 

25.38 

-1.5 

-4.26 


-5.8 


36,000 

24.62 

-2.3 

-4.63 


-6.9 


37,000 

23.92 

-2.6 

-4.99 


-7.6 


38,000 

23.25 

-3.4 

-5.33 


-8.7 


29,000 

22.62 

-3.9 

-5.67 


-9.6 


40,000 

22.02 

-5.1 

-6.00 


-11.1 
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TABLE 5 

INTENSITY CORRECTION INFORMATION 


Percent 

of 

Display 

Width 

Slant 

Range 

(ft) 

Depression 

Angle 

(Degrees) 

ag2 

(dB) 

(dB) 

^eff 

(dB) 

Total 

Error 

(dB) 

0 

16,236 

67.50 


5.75 



4 

17,000 

61.92 


5.65 



10 

18,000 

56.44 

-3.4 

4.40 

-1 .5 

-0.5 

15 

19,000 

52.14 

-2.8 

3.70 

-0.6 

-0.3 

21 

20,000 

48.59 

-2.1 

3.03 

0.0 

-0.9 

26 

21,000 

45.58 

-1.7 

2.39 

0.0 

-0.7 

32 

22,000 

42.99 

-1.3 

1.79 

0.0 

-0.5 

38 

23,000 

40.71 

-1.1 

1.21 

0.0 

-0.1 

43 

24,000 

38.68 

-0.9 

0.65 

0.0 

-0.3 

49 

25,000 

36.87 

-0.6 

0.12 

0.0 

-0.5 

50 

25,236 

36.47 

-0.5 

0.00 

0.0 

-0.5 

54 

26,000 

35.23 

-0.3 

-0.39 

-0.1 

-0.8 

60 

27,000 

37.75 

-0.2 

-0.88 

-0.5 

-1.6 

65 

28,000 

32.39 

-0.1 

-1.35 

-1.1 

-2.6 

71 

29,000 

31.14 

0.0 

-1.81 

-2.0 

-3.8 

76 

30,000 

30.00 

0.0 

-2.25 

-3.0 

-5.3 

82 

31,000 

28.94 

0.0 

-2.68 

-4.5 

-7.2 

86 

32,000 

27.95 

0.0 

-3.09 


-7.6 

93 

33,000 

27.04 

-0.1 

-3.49 


-7.1 

99 

34,000 

26.12 

-0.2 

-3.88 


-8.6 

100 

34,236 

25.98 

-0.2 

-3.97 


-8.7 


35,000 

25.38 

-0.4 

-4.26 


-9.2 


36,000 

24.62 

-0.4 

-4.63 


-9.5 


37,000 

23.92 

-0.5 

-4.99 


-10.0 


38,000 

23.25 

-0.6 

-5.33 


-10.4. 


29,000 

22.62 

-0.7 

-5.67 


-10.9 


40,000 

22.02 

-0.8 

-6.00 


-11.3 
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TABLE 6 

L„,, INTENSITY CORRECTION INFORMATION 

H V 


Percent 

of 

Display 
. Width 

Slant 

Range 

(ft) 

Depression 

Angle 

(Degrees) 

(dB) 

(dB) 

^eff 

(dB) 

Total 

Error 

(dB) 

0 

16,236 

67.50 


5.75 

-10.0 


4 

17,000 

61.92 


5.65 

-10.0 


10 

18,000 

56.44 

-4.9 

4.40 

-10.0 

-10.5 

15 

19,000 

52.14 

-3.8 

3.70 

-9.0 

-9.1 

21 

20,000 

48.59 

-2.9 

3.03 

-8.0 

-7,9 

26 

21,000 

45.58 

-2.4 

2.39 

-7.0 

-7,0 

32 

22,000 

42.99 

-1.8 

1.79 

-6.0 

-6.0 

38 

23,000 

40.71 

-1.4 

1.21 

-5.0 

-5.2 

43 

24,000 

38.68 

-1.1 

0.65 

-4.0 

-4.5 

49 

25,000 

36.87 

-0.8 

0.12 

-2.8 

-3.5 

SO 

25,236 

36.47 

-0.6 

0.00 

-2.5 

-3.1 

54 

26,000 

35.23 

-0.5 

-0.39 

-1.9 

-2.8 

60 

27,000 

37.75 

-0.4 

-0.88 

-1.4 

-2.7 

65 

28,000 

32.39 

-0.2 

-1.35 

-1.0 

-2.6 

71 

29,000 

31.14 

0.0 

-1.81 

-0.8 

-2.6 

76 

30,000 

30.00 

0.0 

-2.25 

-0.7 

-3.0 

82 

31,000 

28.94 

0.0 

-2.68 

-0.5 

-3.2 

86 

32,000 

27.95 

0.1 

-3.09 


-3.7 

93 

33,000 

27.04 

-0.2 

-3.49 


-4.2 

99 

34,000 

26.12 

-0.4 

-3.88 


-4.8 

100 

34,236 

25.98 

-0.6 

-3.97 


-4.9 


35,000 

25.38 

-0,6 

-4.26 


-4.9 


36,000 

24.62 

-0,6 

-4.63 


-5.2 


37,000 

23.92 

-0.7 

-4.99 


-5.7 


38,000 

23.25 

-0.8 

-5.33 


-6.1 


29,000 

22.62 

-0.9 

-5.67 


-6.6 


40,000 

22.02 

-1.0 

-6.00 


-7.0 
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FIGURE 7. DATA-COLLECTION GEOMETRY 
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Radar images with all four intensity corrections would appear 
more uniform over a range of depression angles if the target scene 
were all horizontal with no relief. However, point target intensi- 
ties would be distorted. For the purposes of this analysis, all four 
intensity corrections are considered necessary. 

Figures 8, 9, 10, and 11 show the images which have been 
intensity-corrected (in addition to being SR-to-GR and range-azimuth- 
scale-corrected). The symbol “GSC" is used to show that all these 
corrections have been made. Intensity correction is signified by 
the "C" (for "calibration"). 

It should be pointed out that these data (as well as all other 
image data discussed and presented in this report) are amplitude data 
and not intensity data. If these amplitude data were squared to ob- 
tain intensity data, the resulting dynamic range could not be fully 
represented in digital form because of computer limitations. In 
order to avoid this loss of information, amplitude data were used in 
the digital analysis described in the following chapter. 
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FIGURE 8. 


Xhh eRdc imagery corrected for 

SR-TO-GR, RA.-AZ. SCALE, AND 
INTENSITY. 
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FIGURE 9. Xhv ERDC IMAGERY CORRECTED FOR 
SR-TO-GR, RA.-AZ. SCALE, AND 
INTENSITY. 


1-11 - 
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FIGURE 10. Lhh ERDC IMAGERY CORRECTED FOR 
SR-TO-GR, RA.-AZ. SCALE, AND 
INTENSITY. 
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FIGURE 11. Lhv ERDC IMAGERY CORRECTED FOR 
SR-TO-GR, RA.-AZ. SCALE, AND 
INTENSITY. 
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3 

IMAGE ANALYSIS AND INTERPRETATION 

Analysis of the X-L SAR imagery of the SP Mountain area included 
the visual interpretation of the corrected four-channel SAR data and 
evaluation of two digital image enhancement techniques. Since the 
operating wavelengths were 3 cm (X-band) and 23 cm (L-band), dif- 
ferent roughness scales should be apparent in the different images. 
These comparisons are presented in the following section (3.1). 
Next, image differencing was accomplished to accentuate the differ- 
ential response in the like- and cross-polarized imagery (Section 
3.2). Finally, a method for correction for the ground relief was 
developed and used to correct for distortion caused by vertical 
terrain displacement (Section 3.3). 

3.1 IMAGERY INTERPRETATION 

Interpretation of the SP Mountain SP flow SAR imagery was accom- 
plished with the aid of a Geological Survey map (1962). Two distinct 
surface units are indicated; unit 1 consists of basaltic blocks and 
unit 2 corresponds to sand-filled areas with grassy vegetation cover. 
Some more highly vegetated areas are indicated near the southern fan. 

In general, unit 1 and unit 2 do not reflect much differently on 
the four channels. They are better differentiated at L-band than at 
X-band, perhaps due to the reflectivity of the grass at X-band. The 
northern fan is shown as sandy on its northeast portion, but the 
radar images do not show this. Some of the vegetated areas appear 
as prominent return areas on all four channels, but they do not ap- 
pear significantly different from returns from rocks. 

Detailed comparison of imagery from the various channels points 
out specific features with backscatter signatures that change dra- 
matically with changes in polarization or frequency. Referring to 
Figures 8-11, the following observations can be made: 
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There is a prominent rectangle east of SP Mountain on but 

it's shape is quite different on and it does not show on 
or- 

There is an area of strong return at the base of SP Mountain 
(ESE) on This is fairly prominent in L^y, and barely visi- 

ble on the X-band channels. 


There is a low return strip near the base of SP Mountain (N, E, 
and possibly W) on the L-band channels. This strip is not visible 
on the X-band channels. 

The X-band channels show what may be a road running east and west 
just north of SP Mountain. This does not show on the L-band 
channels . 

Michilbach tank (NE of the flow) shows on all four channels. The 
X^y image does not show the fences both east and west of the tank, 
although they appear on the other three channels. 

There are long lines both north and south of Michilbach tank that 
may be fences. These show brightly at L-band, moderately to weakly 
at Xj^^, and cannot be distinguished at X^y. 


3.2 DIFFERENCING 

One of the ways that two different data channels can be compared 
is by creating a difference image, in which each pixel is the dif- 
ference between the corresponding pixels in the two channels being 
compared. This standard MSS enhancement technique was applied to 
the L-band HH and HV data channels using the ARIES digital image 
processing facility. Further information about the digital computer 
work involved is provided in Appendix B. 

The composite difference image in Figure 12 shows L^^ compo- 
nents in black and L^y in white. Image components that were equal 
in both channels would be expected to cancel and appear grey in the 
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FIGURE 12. DIFFERENCE IMAGERY. 
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composite. Image components that are stronger in one channel would 
be expected to appear in the composite as black or white components, 
depending upon which channel is stronger. Note that processing 
streaks in both channels also appear in the composite. 

Low return areas in both channels appear as grey areas in the 
composite. The low-level image texture is less visible in the com- 
posite image. The differencing operation has evidently caused the 
weak areas to cancel to some extent, leaving a nearly featureless 
grey area. 

The moderate and strong components of both channels are all 
visible but displaced because of poor registration. The rough sur- 
face of the lava flow has texture in both channels and in the com- 
posite. The shorter-range portion of the lava appears to have more 
white (cross-polarized) components than black (parallel). This in- 
cludes the steep near edges of the flow. The far-range portion of 
the flow has more dark (parallel) components. 

The lack of alignment of the two channels is quite noticeable 
throughout the composite image. Strong components of both channels 
can be seen (L^^ dark and light) with a slight displacement 
in range. For example, a single point image within the rectangle 
and another southwest of the rectangle both appear as displaced black 
and white images. The strong linear images northwest of center can 
be seen clearly as displaced pairs also. If these pairs were pro- 
perly registered, they would tend to cancel each other out, rather 
than appearing prominently. A notable exception is a rectangular 
array of components in the southeast portion of the image that 
appears almost exclusively in the cross-polarized channel. It 
appears white in the composite. 

The observed displacement is greater at near range than at far 
range; this non-uniformity may be due to slightly different sweep 
speeds in the two channels. Co-registration of SAR channels is 
therefore probably required, in addition to the slant-to-ground 
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range and aspect ratio correction (see Section 2.3), to produce a 
satisfactory composite image. Identical channels could be used with 
the subtraction routine to test the accuracy of the registration. 


Use of the channel subtraction technique does appear to enhance 
some features of the SP flow and to increase contrast between low 
and moderate to strong returns. Its effectiveness is significantly 
degraded, however, by mi s-registration of the two channels. This 
example points out the need for developing algorithms to co-register 
SAR channels. 

3.3 REMOVAL OF RELIEF DISPLACEMENT 

In a sidelooking radar (SLAR) image, the slant-range representa- 
tion causes two basic geometric distortions. The first is the dis- 
tortion in representing a horizontal terrain segment ("slant-to- 
ground" distortion). The second is the distortion in representing 
vertical terrain segments ("relief displacement"). 

The first kind of distortion is easily removed for flat surfaces, 
such as oceans and midwest agricultural fields. This correction has 
been discussed in Section 2.3.1. 

Correction of the second type of distortion, however, has not 
been accomplished to our knowledge. Relief displacement, illustrated 
in Figure 13, severely distorts SLAR images of mountain and hilly 
terrain; a relative displacement with elevation occurs toward the 
ground track of the aircraft. Such a displacement precludes accurate 
overlays and computer comparison with other imagery such as Landsat, 
causes "radar foreshortening," and, as shown by the relative posi- 
tions of A and B in Figure 13, radar "layover." This latter distor- 
tion occurs when the downward slope of the terrain toward the air- 
craft is greater than the perpendicular to the radar beam. 

Relief displacement can be removed using topographic data which 
are available in digital form for the continental United States, and 
can be obtained from stereo aerial photography elsewhere. 
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radar foreshortening. Note that 
elevation causes relative displacement 
toward the vehicle.) 
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If one commences with the map coordinates (ground range) and the 
flight parameters of the radar, the corresponding slant ranges can 
easily be found. The image intensity at the corresponding slant 
range can then be assigned to the proper map coordinate. 

The procedure for each ground range sample is as follows: 

1. Select a map coordinate (ground range), g, related to the 
known ground track of the vehicle. 

2. Using the known elevation at that coordinate, find the slant 
range s corresponding to g: 

s(g) = (g2 + (h - 

where e is the elevation and h is the altitude of the 
aircraft. 

3. Assign the image intensity found in the corresponding slant 
range s(g) to the map coordinate g. As the slant range is 
discretely sampled, assign the intensity found by nearest 
neighbor, linear interpolation, or a more sophisticated 
interpolation as desired or required. 

The above procedure is presented in its simplest form. Flight 
parameter errors, scaling, and bookkeeping difficulties (such as re- 
lating digital topographic data skewed with SAR data) will provide 
sufficient difficulty in execution. 

Relief displacement in SAR imagery of SP Mountain in northern 
Arizona (contours shown in Figure 14) has been removed. SP Mountain 
is not perfectly rotational ly symmetric, but it is close enough to 
provide a graphic illustration of relief displacement and its re- 
moval. The 3-cm (X-band) cross-polarized (HV) imagery is displayed 
in four manifestations in Figures 2, 15, and 16. 
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FIGURE 14. TOPOGRAPHIC CONTOURS OF SP MOUNTAIN 

(Note deviations from circular 
symmetry and higher elevations on 
west side.) 
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FIGURE 15. SLANT- TO-GROUND-RANGE CONVERSION WITH 
DATUM PLANE AT TOP OF SP MOUNTAIN. 



FIGURE 16. REMOVAL OF RELIEF DISPLACEMENT. 

(Includes slant-to-ground-range 
conversion. ) 
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The slant-range uncorrected imagery shown in Figure 2(a) demon- 
strates slant-range distortion in that both the base of the mountain 
and its crater are foreshortened in the range dimension. Relief 
displacement causes the crater to be imaged unnaturally close to the 
ground track side of the mountain. 

Distortion removal by conventional slant-to-ground range con- 
version is shown in Figure 2(b), where the datum plane is taken at 
the base of the mountain. Here, the base of the mountain is rec- 
tified to approximate its actual shape (see Figure 14), but the 
crater is still displaced in position and unnaturally foreshortened. 

If the datum plane for conventional slant-to-ground range conver- 
sion is taken at the elevation of the crater, the crater's actual 
shape and position are approximated, but the base of the mountain is 
distorted and displaced, as shown in Figure 15. 

Using the procedure described above, both relief displacement 
removal and slant-to-ground range conversion are accomplished. The 
elevation data for SP mountain were approximated by four analytic 
functions for various sections of the scene. As shown in Figure 16, 
the preliminary result is less than ideal, but the actual shape of 
the entire mountain is approximated, with both the shapes of the 
crater and the base of the mountain closely approximated. Additional 
effort is required to refine the topographic approximation and 
develop a more accurate determination of flight parameters. 

Thus, a simple procedure for the removal of relief displacement 
in sidelooking radar imagery has been described and a very prelimi- 
nary test shown. The procedure requires elevation data and knowledge 
of flight altitude and ground track. Fortunately, digital elevation 
data for the continental United States are readily available, and 
can be obtained elsewhere (with more effort and expense) from stereo 
aerial photography. The removal of relief displacement will enable 
accurate overlays of radar imagery upon maps and other types of 
imagery, such as Landsat. 

The computer programs used are discussed in Appendix B. 
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4 

CONCLUSIONS AND RECOMMENDATIONS 


4.1 CONCLUSIONS 

In carrying out the work discussed in this report, ERIM has 
demonstrated the use of digital image processing routines that can 
aid in geological interpretation of SAR imagery. Standard rectifi- 
cation and registration routines have been applied to all four chan- 
nels of SAR SP mountain data and the resulting imagery has been 
manually interpreted. In addition an MSS differencing technique for 
image enhancement has been applied, and an advanced rectification 
•routine for relief displacement has been developed. 

Basic preprocessing routines that should be applied previous to 
interpretation or additional image enhancement include radiometric 
rectification and image registration. The intensity correction of 
the data in the range direction has been used in this program to 
correct for (1) range falloff (R ), (2) the antenna pattern, and 
(3) CRT-photo-recorder variations. Intensity variations in corrected 
imagery correspond to actual variations in radar backscatter. 

Image registration is necessary to compare SAR imagery with maps 
in ground coordinates. These routines are (1) the slant-range-to- 
ground-range correction and (2) the range-azimuth scale correction. 
For the latter, a relatively inexpensive approach was employed which 
permits the modification of the range scale only; this approach gives 
a unity ratio of range and azimuth scales but does not assure exactly 
the same image size. If needed in the future, a more sophisticated 
(and much more expensive) approach is available at ERIM which permits 
changing both range and azimuth scales; if this approach were used, 
exact image registration would be possible. 

Manual interpretation of the corrected imagery showed that the 
L-band channels gave a better representation of the general geologic 
units of the area than did X-band. It was suggested that increased 
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reflectivity of the vegetation at X-band masked rock (or soil) sig- 
nature differences. Some vegetated areas showed strong returns on 
all four channels. Several features showing distinctly different 
signatures in the various polarization or frequency channels were 
also observed. 

In an attempt to enhance observed polarization differences, an 
MSS differencing technique was applied to the L-band channels. This 
technique generates a new image in which each pixel is the difference 
between the corresponding pixels in the two single channels. Al- 
though differencing did enhance the contrast of the scene and make 
some flow features more prominent, accurate registration of the two 
channels proved to be a critical factor in the effectiveness of this 
technique. 

Additional image improvement was obtained with a correction for 
relief displacement. This technique has been developed to remove 
the effects of range layover so that, for example, the crater at the 
top of SP mountain is correctly aligned with its actual geographic 
position. The feasibility of this technique has been demonstrated, 
but additional refinement of the algorithm is required to obtain 
completely satisfactory results. Accurate relief displacement cor- 
rection and registration routines are needed prior to combining SAR 
imagery with maps or Landsat imagery in a digital format. 

4.2 RECOMMENDATIONS 

The results of this analysis have provided a valuable basis from 
which to determine future research projects related to digital image 
analysis. The following recorranendations are direct consequences of 
this study: 
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1. Develop improved registration techniques for SAR data chan- 
nels. The image registration algorithm could be tested by 
using identical data as inputs to the image differencing 
algorithm. If the images are in perfect registration, the 
difference will be zero. Any departure from zero would be a 
measure of lack of registration. 

2. Investigate the application of other MSS image enhancement 
techniques, such as ratioing, to SAR digital data. 

3. Develop an improved relief displacement correction algorithm. 
This will involve refining the accuracy of the topographic 
approximation equations and developing a more accurate 
determination of flight parameters. 

4. Develop a technique to register four-channel SAR digital data 
with Landsat data. 
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DIGITIZATION PARAMETERS 
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DIGITIZATION PARAMETERS 
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Figures A-1 , A-2, A-3, and A-4 show the computer printouts of 
the digitization parameters for tapes IPL 584, 585, 586, and 587. 
The CCTs and these printouts were generated by ERIM's Hybrid Image 
Processing Facility (HIPF). 
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R BCOMP 
«SDSP 

5AR DIGITIZING SYSTEM PARAMETERS 
CURRENT DATE'?220CTe0 

FLIGHT DATE. PASS NUMBER, AUX DATA'i'200CT79 PASS 3 SP MOUNTAIN 
RADAR RANGE SCALE FACTOR. RADAR AZIMUTH SCALE FACTOR 
’’196700,15000 

SIGNAL FILM ASPECT RATIO IS 13.1133 

OPTICAL AZIMUTH DEMAGNIFICATI0N713 . 7 

DESIRED IMAGE RESOLUTION} RANGE. AZIMUTH:'?3.3 

SAMPLES PER RANGE RESOLUTION ELEMENT (NY0UIST=2 ) ’’2 

SAMPLES PER AZIMUTH RESOLUTION ELEMENT (NYQUIST=2)72 

RANGE SCALE FACTOR= 1.5 METERS/PIXEL 

AZIMUTH SCALE FACTOR= 1.5 METERS/PIXEL 

RANGE IMAGE PIXEL SPACING IS 7.62582 MICRONS 

DISSECTOR RELAY LENS MAGNIFICATION'?! 

RANGE SCAN APERTURE IN RADAR SPACE = 3.5406 METER GAUSSIAN 
AZIMUTH IMAGE PIXEL SPACING IS 7.29927 MICRONS 
AZIMUTH IMAGE SLIT UIDTH (MICR0NS)-?5 

AZIMUTH SCAN APERTURE IN RADAR SPACE = 1.0275 METER RECT 
ALONG TRACK START POSITION (KIL0METERS)79.5> 

RANGE OFFSET TO START POINT (KILOMETERS)'?! .3 . i 
NUMBER OF PIXELS PER LINE. NUMBER OF LINES. POINT SKIP 

■?2650.. 6800,0 I 

SWATH WIDTH IS 3.975 KILOMETERS i 

ALONG TRACK LENGTH IS 10.2 KILOMETERS ’ 

START POINT 175 t 

STOP POINT 2825 i 

LASER WAVELENGTH (NAN0METERS?515 i' 

FREQUENCY PLANE WEIGHTING FILTER'?NO 

RANGE FREQUENCY PLANE APERTURE IS 18.9094 MILLIMETERS 

AZIMUTH FREQUENCY PLANE APERTURE IS 1.442 MILLIMETERS I 

DEFLECTION VOLTAGE AT FULL SCAN (3600) IS 2.05623 VOLTS 

DEFLECTION VOLTAGE AT STOP POINT IS 1.29738 VOLTS j 

LASER POWER (MW). OBJECTIVE FL (MM). MULTIPLIER VOLTAGE. GAIN I 

7400.8,1800.5 

CONVERSION MODE (LINEAR OR SQRT)7SaRT 

INTEGRATION TIME (MICROSECONDS ) ?990 

SIGNAL FILM RUN LENGTH IS 680 MILLIMETERS | 

SIGNAL FILM DRIVE FREQUENCY IS 388.683 HERTZ 
FILM DRIVE SIGNAL PERIOD IS 2572.79 MICROSECONDS 
TIME PER LINE IS. 2. 76086 SECONDS I 

TUN TIME IS 312 MINUTES 

I'EET OF RECORDING TAPE REQUIRED = 2302 I 

1792.93 4553.31 

P f^S3 7^ 

DO . 

i ABEL7IPL 584 FILE 1 FLIGHT 200CT79 X-BAND HH SP MOUNTAIN 2650X6800 1.5M/PIX SORT 22DCT80 
' 'f'ES'?6800 
■ ■ UO 
' I'JF 

: U TOTAL 7 PAD LINES 


FIGURE A-1 . 


HIPF PARAMETER PRINTOUT FOR CCT NO. 


I PL-584. 
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R BCOMP 
»SDSP 

SAR DIGITIZING SYSTEM PARAMETERS 
CURRENT DATE72Z0CT80 

FLIGHT DATE. PASS NUMBER. AUX DATAT>200CT79 PASS 3 SP MOUNTAIN 
radar range scale factor, RADAR AZIMUTH SCALE FACTOR 
T195200. 15000 

SIGNAL FILM ASPECT RATIO IS 13.0133 I 

OPTICAL AZIMUTH DEMAGNIFICATI0N7 13 . 7 I 

desired image RESOLUTION! RANGE. AZIMUTH! 'i’3 . 3 
SAMPLES PER RANGE RESOLUTION ELEMENT ( NYQUIST=2> fZ 
SAMPLES PER AZIMUTH RESOLUTION ELEMENT < NY0UIST=2) -»2 
range scale FACT0R= 1.5 METERS/PIXEL 
AZIMUTH SCALE FACTOR= 1.5 METERS/PIXEL ' 

range image PIXEL SPACING IS 7.60442 MICRONS 
DISSECTOR RELAY LENS MAGNIFICATIONTl 

range SCAN APERTURE IN RADAR SPACE = 3.5136 METER GAUSSIAN 
azimuth image pixel SPACING IS 7.29927 MICRONS 

AZIMUTH IMAGE SLIT WIDTH ( MICRONS ) ’S j 

azimuth scan APERTURE IN RADAR SPACE = 1.0275 METER RECT I 

ALONG TRACK START POSITION (KILOMETERS )'!’9 . 5 | 

range offset TO START POINT < KILOMETERS IT 1 . 3 i 

number of pixels per line, number OF LINES, POINT SKIP 1 

■'2650,6800.0 | 

SWATH WIDTH IS 3.975 KILOMETERS j 

along track LENGTH IS 10.2 KILOMETERS I 

START POINT 175 
STOP POINT 2825 

LASER WAVELENGTH <NAN0METERSZ5I5 
FREOUENCY PLANE WEIGHTING FILTERTNO 

range frequency plane APERTURE IS 18.7652 MILLIMETERS > 

AZIMUTH FREQUENCY PLANE APERTURE IS 1.442 MILLIMETERS : 

deflection VOLTAGE AT FULL SCAN <3600) IS 2.07203 VOLTS ' 

DEFLECTION VOLTAGE AT STOP POINT IS 1.30735 VOLTS ' 

LASER POWER (MW). OBJECTIVE FL (MM), MULTIPLIER VOLTAGE, GAIN i 

^400.8.1800.5 1 

CONVERSION MODE (LINEAR OR SQRT)''SQRT 
INTEGRATION TIME (MICROSECONDS ) ''990 
SIGNAL FILM RUN LENGTH IS 680 MILLIMETERS 
SIGNAL FILM DRIVE FREQUENCY IS 388.683 HERTZ 
FILM DRIVE SIGNAL PERIOD IS 2572.79 MICROSECONDS 
TIME PER LINE IS 2.76096 SECONDS 
RUN TIME IS 312 HINUIES 
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FIGURE A-2. HIPF PARAMETER PRINTOUT FOR CCT NO. IPL-585. 
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FIGURE A-3. HIPF PARAMETER PRINTOUT FOR CCT NO. IPL-586. 
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B.1 INTRODUCTION 

The following image data for four channels in the slant-range 
plane were provided: 

X Band HH Polarization - Tape IPL 584 

X Band HV Polarization - Tape IPL 585 

L Band HV Polarization - Tape IPL 586 

L Band HH Polarization - Tape IPL 587 

These data were supplied in byte form and the image sizes (no. of 
records x no. of elem./rec.) are given in Table B-1 . 

TABLE B-1 
IMAGE SIZES 


Image 

Size 

IPL 584 

6800 X 2650 

IPL 585 

6800 X 3000 

IPL 586 

6800 X 3000 

IPL 587 

6800 X 3000 


B.2 IMAGE RE-SCALING 
B.2.1 FIRST APPROACH 

The first part of this task requires re-scaling of the images, 
both in range and azimuth. 

The software available for this operation in the ARIES facility 

is a re-sampling algorithm under UIC[141,21] and task name DSMPRD. 

★ 

This program requires complex-word data as input and a program 
was written to reformat the image data. 

*UICL144,21J Task RTOC. 
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c 

The size of the newly-formatted data is a huge 81.6 x 10 
bytes, and even in a large disk, it is difficult to find space for a 
contiguous file of this size. Consequently, the complex-valued for- 
matted data was generated in three sections. 

Input data buffer size in the present POP 11/45 system is re- 
stricted to a maximum 2048 complex words in the re-sampling program 
(DSMPRD). Accordingly, each record must be processed in two 
sections. 

In conclusion, this procedure requires six program runs for each 
channel . 


B.2.2 SECOND APPROACH 

The procedure above is very time-consuming and it does not seem 
to be cost-effective. The azimuth scale change requirement was 
therefore dropped and the range scale factor was adjusted to maintain 
the same aspect ratio. 

The scale change (in the range direction) can now be accomplished 
with a program designed to convert data from the slant-range plane 
to the ground-range plane with a range scale change capability. This 
software is implemented in the ARIES facility in UIC[150,033] under 
the task name STG. 

The histograms corresponding to the four tapes generated are 
shown in Figures B-1 , B-2, B-3, and B-4. 

B.3 INTENSITY CALIBRATION 

The purpose of the calibration is to affect the radar data by 
the antenna polar pattern and the recorder transfer function multi- 
plicative corrections in order to obtain an accurate measurement of 
the terrain reflectivity at a particular wavelength and polarization. 
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FIGURE B-2. HISTOGRAM FOR TAPE I PL-585 GS (X, 
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FIGURE B-3. HISTOGRAM FOR TAPE I PL-587 GS (L 
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FIGURE B-4. HISTOGRAM FOR TAPE I PL- 586 GS (L 


RADAR AND OPTICS DIVISION 



This calibration can be implemented in the slant-range plane 
(preferably) or in the ground-range plane. This latter implementa- 
tion should be affected by any possible scale changes previously 
applied. 

In the ARIES facility, a software package* was implemented 
that reads a small number of samples of the calibration vector, fits 
a polynomial function to these data (least-squares criterion), and 
evaluates the calibration vector components for all range locations. 

The histograms corresponding to the four tapes generated are 
shown in Figures B-5, B-6, B-7, and B-8. 

B.4 IMAGE DIFFERENCE 

The image difference between the L-band HH and HV channels was 
obtained. This difference was obtained pixel-by-pixel according to 
the expression 

DIF = IMl - IM2(E1/E2) + BIAS, 

where El and E2 are the average values of image 1 and image 2, re- 
spectively. A bias was used to allow an accurate representation of 
the bipolar nature of the image. 

Given the slightly different scales of image 1 and image 2, the 
difference obtained is not totally accurate and thus is enhanced 
along edge-type features. 

Two software packages were implemented in the ARIES facility to 
obtain the image difference. A slow Fortran task [144,21]INDIF.TSK 
and a fast AP120B based [144,21]1MDIF.APP were used. Both read two 
disk-based images in byte form and generated to disk a byte-form 
image difference. See Figure B-9. 

Figure B-10 shows the histogram for the difference image data. 

*UIC[144,21J Tasks Read and Inter. 
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FIGURE B-5. HISTOGRAM FOR TAPE I PL-584 GSC (X, 


RADAR AND OPTICS DIVISION 
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FIGURE B-6. HISTOGRAM FOR TAPE IPL-585 GSC (X, 


2TRi 


RADAR AND OPTICS DIVISION 
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FIGURE B-7. HISTOGRAM FOR TAPE IPL-587 GSC (L 


RADAR AND OPTICS DIVISION 
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FIGURE B-8. HISTOGRAM FOR TAPE I PL-586 GSC (L 
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FIGURE B-10. HISTOGRAM FOR IMAGE DIFFERENCE DATA. 


RADAR AND OPTICS DIVISION 


b.5 correction of relief displacement 

The implementation of the algorithm to remove relief displacement 
from the SP Mountain Data was carried out and its feasibility was 
shown. The elevation model utilized is, necessarily, inaccurate and 
this is the reason for the less-than-perfect reconstruction. 

The ARIES software developed and used for this purpose is the 
following: 

1. [144,21]C0NT0UR.TSK - A Fortran program modeling the eleva- 
tion map of SP Mountain. See Figure B-11. 

2. [144,21 ]TNSTOG.TSK - A Fortran program that reads disk-based 
slant-range image data and a disk-based elevation map and 
then reconstructs a ground-range map with elevation cor- 
rection. This particular program does not achieve a unity 
aspect ratio. This is due to the fact that (normally) the 
sample spacing is the same in range and azimuth for the 
slant-range data. Upon SR-to-GR conversion, the ground-range 
data sample spacing changes and, to a first-order approxi- 
mation, it is assumed constant along the range of interest. 
A second display program would be required to achieve unity 
aspect ratio. Given the initial stage of this idea, an 
approximate unity aspect ratio was obtained using a simple 
technique. See Figure B-12. 

Figures B-13 and B-14 show the histograms for the relief dis- 
placement correction data. Figure B-13 corresponds to the SR-to-GR 
correction while Figure B-14 corresponds to the complete correction. 
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FIGURE B-n. PROGRAM CONTOUR FOR RELIEF DISPLACEMENT CORRECTION. 
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B-16 


RADAR AND OPTICS DIVISION 



rv TC 5u 
C 
C 
C 

C FI.'IC rUilV^TIOr’S T’4 CFi^TFAL CKATEF 

C 

c 

10 = 

FLFV = u-\/h-3UPT(r-**2-nST»*'? ) 

F ( I ) =M'-T I ELtV ) 

0 TO 5 u 

20 IF (i;i',T.ijr.hH(-)GC TO ■+0 

2U EL EU=A-\/C J.ST 

Em=f'I'T(CLEV) 

GO TC so 

“0 F.LEV=h2C0. 

E(I)=MMT(ELaM 
5u COMIkUc 

C 

c 

JCUTru 

CALL '.<RITE(CTFL1.ERP(3) .t.HE.uCLiT) 
CALL rO-AiT(CTPLl.e;RR<6) ) 

c 

c 

lOno COi'iTlfOc 

c 

c 

c 

CALL CETfiCh(CTRLl.IER«) 

C 

C 

C 

EI.l 

c 

c 

c 
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PROGRAM CONTOUR FOR RELIEF DISPLACEMENT CORRECTION. 
(Concluded) 
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FIGURE B-12. PROGRAM TNSTOGE FOR RELIEF DISPLACEMENT CORRECTION. 
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FIGURE B-12. PROGRAM TNSTOGE FOR RELIEF DISPLACEMENT CORRECTION. 

(Continued) 
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FIGURE B-12. PROGRAM TNSTOGE FOR RELIEF DISPLACEMENT CORRECTION. 

(Concluded) 
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FIGURE B-13. HISTOGRAM FOR SR-TO-GR CORRECTION DATA. 
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FIGURE B-14. HISTOGRAM FOR SR-TO-GR AND RELIEF DISPLACEMENT CORRECTION DATA. 


